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Abstract 

The emission from blazars is known to be variable at 
all wavelengths. The flux variability is often accom- 
panied by spectral changes. Spectral energy distri- 
bution (SED) changes must be associated with changes 
in the spectra of emitting electrons and/or the phys- 
ical parameters of the jet. Meaningful modeling of 
blazar broadband spectra is required to understand 
the extreme conditions within the emission region. 
Not only is the broadband SED crucial, but also in- 
formation about its variability is needed to under- 
stand how the highest states of emission occur and 
how they differ from the low states. This may help 
in discriminating between models. Here we present 
the results of our SED modeling of the blazar S5 
0716-1-714 during various phases of its activity. The 
SEDs are classified into different bins depending on 
the optical brightness state of the source. 
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1. Introduction 

S5 0716-1-714 is a bright, high declination BL Lac 
object at a redshift, z = 0.31 ±0.08 (Nilsson et al. 
2008). This source has been extensively studied across 
the whole electromagnetic spectrum and exhibits strong 
variability on a wide range of timescales, ranging 
from minutes to years (e.g., Wagner et al. 1990; Heidt 
& Wagner 1996; Villata et al. 2000; Raiteri et al. 
2003; Montagni et al. 2006; Ostorero et al. 2006; 
Gupta et al. 2008a, b, 2009 and references therein). 
Nearly periodic oscillations of '--'15 minutes in the 



optical R band were detected in this source (Rani et 
al. 2010). The optical duty cycle of S5 0716-1-714 
is nearly unity, indicating that the source is always 
in an active state in the visible (Wagner & Witzel 
1995). This blazar was recently shown to be a strong 
source in the high energy gamma-ray band by Fermi- 
LAT (Abdo et al. 2009). 



2. Multi-frequency Data 

We carried out this study of the SEDs of the BL 
Lac object S5 0716+714 using high quality multi- 
frequency data in the literature. The data cover ra- 
dio to optical energy bands. The optical U, B, V, R 
and I fluxes were collected from Raiteri et al. (2003), 
Montagni et al. (2006), Villata et al. (2000), Os- 
torero et al. (2006) and Gu et al. (2006). The data at 
radio frequencies of 22 and 37 GHz were taken from 
Salonen et al. (1987), Terasranta et al. (1992, 1998, 
2004, 2005). The UMRAO^ (University of Michi- 
gan Radio Astronomy Observatory) data at 4.8, 8 
and 14.5 GHz frequencies was provided by Margo 
AUer. The data at optical frequencies spans over a 
period of 1994 November through 2004 April while 
the radio data cover a time period of 1991 September 
to 2005 May (Fig. 1). 

We model the SEDs of S5 0716-^714 using the 
above radio to optical frequency data. The long term 
data allow us to obtain six different mean SEDs of 
the source corresponding to different phases of its 
activity (Fig. 2). These six different SEDs are char- 
acterized by different optical outburst phases and the 
radio data is averaged over the same time periods. 
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Figure 1. The multi-frequency radio and optical data of 
the BL Lac object S5 0716+714. 

3. Analysis and Results 

3.1 SED Modeling 

We use a homogeneous synchrotron self-compton 
model (SSC) with a broken power law (BPL) to fit 
the lower energy part of observed spectra which are 
characterized by synchrotron spectra. The SED model 
fitting is achieved by using a SED code available on- 
line^. The best fit model was obtained by varying 
the parameters of a numerical SSC code (Tramacere 
2007; Tramacere et al. 2009). The model assumes 
that radiation is produced within a single zone of 
jet (~ radius R), moving relativistically at a small 
angle to the line of sight of the observer. The ob- 
served radiation is amplified via the boosting factor 
6=[T{1- f3 cose)]-K 

Below and above the peaks, the spectrum can be 
approximated with power law profiles with indices 
ai and Q2, respectively. The power law spectra in 
AGNs are naturally produced if the emitting elec- 
tron follow a power-law distribution in energy. We 
approximate this behaviour with a broken power law 



Figure 2. The optical R passband flux of S5 0716+714 
divided into six different bins on the basis of its optical 
activity. Binl, Bin4 and Bin6 represent the optical out- 
burst phases while Bin2, Bin3 andBinS represent dimmer 
phases of the source. 

(BPL) with indices ni and n2, respectively below 
and above the break energy TftmeC^ : 



7\Tp"2-rti-p-n2 
b 



ifr<r6 
ifr>rfe 



(1) 



http://tools.asdc.asi.it/ 



With these approximations, we can completely 
specify the model with the following parameters: maj 
netic field intensity (B), size of emission region (R), 
Doppler boosting factor (6), power-law indices (ni, 
n2), number density of emitting electrons (N) and 
Lorentz factor of the electrons at the break energy 
(Tb). The values of all these fitted parameters for 
different SEDs are listed in Table 1. 

4. Discussion and Conclusions 

4.1 Limitations to the model 

Although, as seen in Fig. 3, we were able to achieve 
reasonably good fits of the synchrotron emission for 
all six SEDs of the source, one should bear in mind 
that the one-zone BPL model is over-simplified in 
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Figure 3. The SSC modeled SEDs of the source corresponding to six different Bins. 



accounting for the radio-optical blazar emission. The 
applicability of a single-zone emitting region has been 
questioned by a number of authors (e.g., Vittorini et 
al. 2009, Raiteri et al. 2010). They showed that the 
BL Lac SED can be more successfully modeled with 
two synchrotron components (two different emitting 
populations). 

Furthermore, unfortunately, we do not have syn- 
chrotron peak measurements for the the source, which 
would significantly help to constrain the model. An- 
other limitation to the model is that we do not cor- 
rect our optical measurements for any possible host 



galaxy contribution. Some objects, such as BL Lac 
itself, have a significant contribution from starlight 
to the optical bands, which will modify the calcu- 
lated synchrotron emission in this region, especially 
during the low states. We also do not make any 
corrections for Galaxy or internal absorption, which 
may again slightly affect the optical-UV part of the 
spectra. Last, but not least, we stress that our data is 
not strictly simultaneous but has been averaged over 
a period of months. As blazars are highly variable 
over timescales of a day or less, the time differences 
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Table 1. Fitted parameters of SED model 
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Binl 


16.8 


0.4 


12.8 
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0.2 
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0.4 
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3.20 
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-9.69 


Bm2 


16.8 


0.4 
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0.2 
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0.6 
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-9.99 


Bin6 


16.8 


0.4 


16.0 


25 


0.2 


6 


0.2 
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R : Size of emitting region 

B : Magnetic field 

& : Doppler boosting factor 

N : Number density of emitting electrons 

TmirK ^max '■ Minimum and maximum values of Lorentz factor 
Ff, : Lorentz factor corresponding to break energy of electrons 
m, 712 : Power law indices 
'^Peak '■ Synchrotron peak frequency 



and averaging might have compromised somewhat 
the modehng. 



4.2 Spectral Energy Distribution variation 

To attempt a study of how the physical parameters 
related to emission region and synchrotron emission 
are changed when the BL Lac S5 0716-1-714 goes 
through various phases of its activity we used the 
simplest approach by fitting a single-zone SSC model. 
The observed results can be summarized as : 

1 . No change between the R and B bands in the mod- 
eled SED are seen during different phases of activity. 

2. The Doppler boosting factor S is higher during the 
optically bright states of the source compared to the 
dimmer phases of activity. 

3. The number density (N) of electrons emitting 
synchrotron photons is larger when the source is in 
lower states. 

4. The synchrotron peak frequency {vpeak) and peak 
intensity {vFypeak) are comparatively higher during 
the optical outburst phases of the BL Lac object. 

This research was supported by CSIR Foreign 
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